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ABSTRACT 

We present the results of our first year of quasar search in the ongoing ESO public Kilo 
Degree Survey (KiDS) and VISTA Kilo-Degree Infrared Galaxy (VIKING) surveys. 
These surveys are among the deeper wide-held surveys that can be used to uncover 
large numbers of z ^ 6 quasars. This allows us to probe a more common population of 
z ^ 6 quasars that is fainter than the well-studied quasars from the main Sloan Digital 
Sky Survey. From this hrst set of combined survey catalogues covering ^250 deg^ we 
selected point sources down to Zab = 22 that had a very red i — Z {i — Z > 2.2) 
colour. After follow-up imaging and spectroscopy, we discovered four new quasars in 
the redshift range 5.8 < z < 6.0. The absolute magnitudes at a rest-frame wavelength 
of 1450 A are between —26.6 < M1450 < —24.4, conhrming that we can hnd quasars 
fainter than M*, which at z = 6 has been estimated to be between M* = —25.1 and 
M* = —27.6. The discovery of four quasars in 250 deg^ of survey data is consistent 
with predictions based on the z ^ 6 quasar luminosity function. We discuss various 
ways to push the candidate selection to fainter magnitudes and we expect to hnd 
about 30 new quasars down to an absolute magnitude of M1450 = —24. Studying this 
homogeneously selected faint quasar population will be important to gain insight into 
the onset of the co-evolution of the black holes and their stellar hosts. 
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1 INTRODUCTION 

Supermassive black holes (SMBHs) are inferred to reside at 
the nuclei of most (if not all) massive galaxies. In present 
day galaxies, tight correlations are observed between the 
black hole mass and global host parameters over at least 
three orders of magnitude in black hole mass (see e.g 


* Based on observations collected at the European Southern 
Observatory, Chile, programmes 177.A-3016, 179.A-2004, 089.A- 
0290, 089.A-0596, 090.A-0383, 090.A-0642 and 091.A-0421. 
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iKormendv fc h 3|2013| . for a recent review). Going back in 
time, the global history of star formation and of radia¬ 
tive mass accretion on to black holes display similar be- 
haviour as a function of time, peaking at redshifts ^2 — 3 
jHeckman et al.ll2004l : for a review see lMadau &: DickinsonI 
I 2 OI 4 I and references therein). Determining the physical 
processes that govern the symbiotic relationship between 
SMBHs and their host galaxies is an outstanding challenge 
in astrophysics. Studying the SMBHs and stellar hosts of 
high-redshift (z > 5.7) quasars is one of the most di¬ 
rect ways to obtain insight on the onset of this relation 
in the first Gyr after the big bang. They provide funda- 
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mental constraints on the formation and growth of the first 
SMBHs, on early star formation and on chemical enrichment 
in g alaxy hosts and the surronnding intergalactic matter 
(e.g. Fan et al. 20061: Jiang et ^120071 : iDe Rosa et al.llioil] : 
ICarilli fc Walteij|20lk 

Thanks to large imaging snrveys of moderate depth, 
it has become feasible to identify the rare bright qnasars 
photometrically at high redshifts. As an illnstration, a In- 
minous quasar at redshift 6 can be detected photometri¬ 
cally in the 2 band within two minutes on a 2.5 m tele¬ 
scope. In the past decade abont 70 bright qnasars at 2 > 5.7 
have been discovered (e.g. Fan et al.||2006l: Jiang et al.ll2009l: 
WilIottet_^ 201I)b||_ MbrtJock et al.lboil : Bahados et al.l 


20141 : IVenemans et alJ l2015ll using a variety of different 


multiwavelengt h surveys such as th e Sloan Digital Sky 
Survey (SPSS: [Abazaiian et al.ll2009l~l. the Canada-Franc e 
High -2 Quasar Survey (CFHQS: Iwillott et al.l l2010lj) . 
the UK I nfrared Telescope In frared Deep Sky Survey 
(UKIDSS; [Lawrence et al.l l2007f l and Panoramic Survey 
Telescope fc Rapid Response System 1 ( Pan-STARRSl; 
iMorganison et al.ll2012l : iBanados et al.llioi^ L 

These quasars appear to have SMBH masses equiva¬ 
lent to today’s most massive galaxies with high-accretion 
rates, close to Eddington. The masses have been inferred 
from quasar emission line width and continuum luminos¬ 
ity extrapolating empirical relationships established in low 
redshift active galactic nuclei (AGN) through multi - epoch 
reverberation mapping. For example, IPe Rosa et al.l (1201 il l 
obtain for a population of 2 ~ 6 quasars a mean black hole 
mass of log(MBH/MQ) « 9.06, typical bolometric luminosi¬ 
ties > lO’^^ erg s“^ and Eddington ratios Lboi/LEdd~ 0.6. 
The growth time for such black holes and Eddington ratios 
if starting from 10^ Mq seeds is comparable to the age of 
the Universe at 2 ~ 6 (~0.9 Gyr). The quasar emission lines 
also suggest well-developed metallicities. The relative abun¬ 
dance of Fe and a elements, inferred from Fell and MgII 
emis sion lines, appears similar to quasars at redshift 4 
(e.g. Ijiang et al.ll2007l : iDe Rosa et ahllioill l. 

With the advent of ALMA it is now possible to obtain 
constraints on host dynamical masses and on star forma¬ 
tion by observing the cool molecular gas of the hosts. The 
quasar hosts appear smaller or at most roughly consistent 
with the p resent-day relation between host mass and SMBH 
mass (e.g. 2004 20091 : IWang et al.ll201ol . l201,jl : 

IWillott. Bergeron fc Omontll201,^ L These dynamical masses 
are based on line kinematics. Direct detection of the stel¬ 
lar hosts requires higher spatial resolution and has to await 
the James Webb Space Telescope and the next generation 
groun d-based 30m-class telescopes (see e.g. iMechtlev et al.l 
I 2 OI 2 I L The star formation rates measured in the quasar host 
galaxies are tens to hundred solar masses per year (e.g. 
Wang et al.ll200^. l2013l : IWalter et ahllioo^ : I Venemans et al.l 


Wang ( 

2012h . 


The results listed above are mostly based on the bright¬ 
est (‘tip of the iceberg’) of the quasar population. The 
quasars have been selected close to the detection limit of 
flux-limited surveys and represent the bright end of a more 
common fainter population. They were detectable thanks to 
the combination of high SMBH mass with high accretion 
rates. T he few notable excepti ons are the M^t^ SMBH s 
found bv lWillott et al.l (l 2010 al l and iKashikawa et al.l 1I2OI5II . 

Given the steepness of the quasar luminosity function 


(see e.g. Fig.Hl IWillott et '^l2010bl : lKashikawa et ai1l2015l l. 
a 2.5 mag (a factor 10 in continuum flux density) deeper 
photometric selection would open the door to a quasar pop¬ 
ulation which is expected to be an order of magnitude more 
abundant. For a given accretion rate (black hole mass) this 
population would probe active galaxies with an order of 
magnitude lower black hole mass (accretion rate). 


The abundance of SMBHs at lower masses and accre¬ 
tion rates is a valuable constraint for the various alternative 
theories for the formation and early evolution of SMBHs 
and their hosts. These include the remnants of Population 
HI stars, yielding black hole seeds of the order of 10^ Mq, di¬ 
rect collapse of primordial dense gas in protogalaxies leading 
to masses of the order of Mq and run-away dynami¬ 

cal instabilities in dense stellar systems which could result 
in i nitial black hol e masses of the order of 10^“® Mq (see 
e.g. IVolonterill201^ . for a review). The theories predict aver¬ 
age mass densities for SMBHs seeds at 2 S> 6 that differ by 
up to an order of magnitude. Hydrodynamical simulations 
are being used to explore in detail these different formation 
channels. They are now able to predict the relative abun¬ 
dance of bright and faint quasars (e.g. Costa et al. I I 2 OI 4 
iKatz. Siiacki fc Haehnei3l2015l : iLatif et al.ll20ld ) . Thus es¬ 
tablishing the quasar demography (accretion rates, star for¬ 
mation rates, host masses) for 10® Mq SMBHs provides 
a direct discriminant on the simulated formation scenarios. 


The demography of quasars with -^10® Mq SMBH 
masses is also important for a fair comparison between the 
SMBH-host galaxy relations at high redshift and today. For 
today’s SMBH-host relations, selecting black holes which 
are 10 times lower in mass corresponds to probing stellar 
hosts which are 10-100 times more abundant, are -^2 mag 
fainter and have ~35% lower stellar velocity dispersions. At 
black hole masses of 10® Mq we are also not hampered any¬ 
more by the potentially severe selection bias in comparing 
AGN-selected samples at high redshift to host-selected sam¬ 
ples at low redshift. It could be that current high-redshift 
samples are overwhelmed by exceptionally massive 10® Mq 
BHs in the more common smaller hosts compared to the 
muc h rarer massive galaxies with equally m assive black holes 
(e.g. Iwillott et aLll^05l : iLauer et al.ll2007l L 


For these reasons we are building up a homogeneous 
sample of faint quasars at 2 ~ 6 by combining ESO’s pub¬ 
lic surveys the Kilo Degree Survey (KiDS) and the VISTA 
Kilo-Degree In frared Galaxy (VIKING) survey which cove r 
the same area (lArnaboldi et al.]|2007l : fPe Jong et aI.|[201^ L 
The nine-band u through Kg photometry from the combined 
surveys goes up to ~2 mag deeper than SDSS, UKIDSS and 
Pan-STARRSl. Here we report on the high -2 quasar har¬ 
vest from the first year of operations: four quasars in range 
5.8 < 2 < 6.0. The results suggest we can expect to build 
up a homogeneous sample of the order of 30 faint quasars 
at 2 6. 


This paper is organized as follows. In Section [21 we 
present the survey data from KiDS and VIKING we made 
use of and their calibration. The selection of quasar candi¬ 
dates is detailed in Section EH followed by a description of 
the imaging and spectroscopic follow-up in Sections 13.21 & 
El The newly discovered quasars are presented in Section |4] 
and we conclude with a discussion and outlook in Section [S] 

In this paper, we adopt a cosmology with Dm = 0.28, 
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SIa = 0.72 and Ho = 70 kms ^ Mpc ^ (iKomatsu et al.l 
[20 i 3). All magnitudes are given in the AB system. 


2 KIDS AND VIKING SURVEYS AND DATA 


KiDS and VIKING a re sister surveys that cover the same 
1500 deg^ on the sky llArnaboldi et al.lboOTl : IPe Jong et al.l 
l201,'j l. The survey areas are divided over two main strips lo¬ 
cated in the Northern and Southern Galactic Cap (NGP 
and SGP). The two strips, roughly equal in area, span 
-^10° in declination centred on (a, 5)=(12.5h,0°) and 
(a,(5) = (0.8h,—30°). A smaller strip is centred on the 
GAMA09 region with (tr, (5) = (9h,0°). KiDS covers the opti¬ 
cal Sloan u, a, r and i ban ds using the OmegaCAM wide- 
held imager I^Kuiikenll2011^ on the VLT Survey Telescope 
(VST) at ESO’s Paranal Observatory. VIKING covers the 
near-infrared Z (hereafter Zy), Y, J , H and Ka using the 
VISTA InfraRed CAMe ra fVIRCAMjPalton et al.ll200dl at 
the V ISTA telescope (lEmerson. McPherson fc Sutherland! 
I 2 OO 6 I I also at Paranal. VIKING and KiDS started regular 
survey operations on 2010 February 15 and 2011 October 
15, respectively. The exposure times and depth of the hlters 
used in these two surveys are summarized in Table [T] 

In this paper we focus on a set of 267 KiDS i band 
slightly overlapping pointings, each with an -^l.lSdeg^ 
held of view. They were observed between the nights of 
2011 September 2 and 2012 June 7, mostly during bright 
Moon. Due to slightly different areal scheduling of KiDS 
and VIKING, the total area with i, Zy and Y band cov¬ 
erage was 254.3 deg^ with 96.1 deg^ covered in the NGP, 
93.2 deg^ in the SGP and 65.0 deg^ was covered in GAMA09. 
The search for high redshift quasars in this early phase of 
KiDS operations was approached as a pilot project for the 
KiDS QSO team. Its aim was to verify the observational 
quality (in terms of photometry, astrometry, artefacts, data 
foibles) for quasar candidate selection and the turn-around 
time between observation and hnal conhrmation via follow¬ 
up spectroscopy. We created preliminary i band calibrated 
imaging and catalogues using the generic wide-held imag- 
ing pipeline in the As tro- WlSlfl survey handling astern 
(iMcFarland et ^l2013| l. as operated at OmegaCEhjj. The 
processing steps applied to the KiDS i band data involved: 


• removing instrumental hngerprint: de-biasing, hatheld- 
ing, preliminary illumination correction and de-fringing, 

• primitive hagging/masking for pixels affected by satu¬ 
ration, cosmic rays and satellite tracks, 

• preliminary photometric calibration (nightly zero- 
points) and astrometric calibration (per-chip in a single ex¬ 
posure) , 

• swarping and conditions of the hve dithers. 

The sou rce extraction was per formed using SEX- 
TRACTOR llBertin fc Arnout3ll996li conhgured for point 
source photometry. The KiDS pointings have overlap 
with the hrst public release of the Kilo-Degree Survey 


^ Astronomical Wide-field Imaging System for Europe, 
http: //www.astro-wise.org 

^ OmegaCEN expertise centre for wide-field imaging, 
http://www.astro.rug.nl/~omegacen 


Table 1. KiDS and VIKING filters, observing time, nominal 
depth and typical seeing conditions. 


Filter 

Ac 

(A) 

Exposure time 

(s) 

Mag limit 
(AB 5 <t 2 arcsec) 

Seeing 

(arcsec) 

u 

3550 

1000 

24.2 

1.0 

9 

4775 

900 

25.1 

0.8 

r 

6230 

1800 

24.9 

0.7 

i 

7630 

1200 

23.5 

0.8 

Zy 

8770 

480 

22.7 

1.0 

Y 

10200 

400 

22.0 

1.0 

J 

12520 

400 

21.8 

0.9 

H 

16450 

300 

21.1 

1.0 

Ks 

21470 

480 

21.2 

0.9 


ijPe Jong et al.ir2015h . However, the calibration used here 
predates the calibration of the public release. 

For VIKING data we used the images an d catalogues 
processed with the VISTA Data Flo w System (iLewis et al.l 
I 2 OO 5 I : iLewis. Irwin fc Bunclarkl2010lf by the Cambridge Sur¬ 
vey Unit (CASU). Catalogues were also retrieved from 
the VISTA Science Archive of the Wide-field Astronomy 
Unit. At the time of our analysis, we could make use of 
VIKING data produced by the vl.1/1.2 CASU-VIKING 
pipeline (’CASU VERS=1.1/1.2’). We refer the reader to 
IVenemans et al.l ll2013h for a description of VIKING and the 
data products used for the selection of quasar candidates. 


3 PHOTOMETRIC SELECTION FOR 

5.8<Z<6.5 QUASAR CANDIDATES AND 
FOLLOW-UP OBSERVATIONS 

3.1 Candidate selection 

The selection of candidate quasars at z > 5.8 was done as 
follows. First, the KiDS i band catalogues were matched 
with the VIKING ZyYJHKs catalogues. VIKING sources 
that had no match in the KiDS catalogues within a 2 arc- 
sec radius were marked as undetected in the i band and a 
3a limiting magnitude appropriate for that corresponding 
KiDS image was set. The median 3a limiting magnitude of 
the KiDS data we analysed was 24.0. Redshift ~ 6 quasar 
candidates were selected in the Zy band with a signal-to- 
noise ratio S/N > 7 and were required to be point sources. 
Point sources were ide ntified as objects with a galaxy proba¬ 
bility Pgai < 0.95 f see I Venemans et al.ll2013l for a discussion 
on the efficiency to select point source with Pgai < 0.95). 
To avoid spurious sources with only a detection in the Zy 
band, we further required that an object was also present 
in the Y band catalogue. We subsequently applied the fol¬ 
lowing colour criteria to select quasars in the redshift range 
5.8 < 2 < 6.4: 

i-Zy> 2.2, Zy - Y < 1.0 and -0.5 <Y - J < 0.5. 

These colou r criteria are base d on th e work by e.g. lFan et al.l 
ll200lll and IVenemans et al.l (l2007l i and are illustrated in 
Fig. [ 1 ] In short, quasars at z > 5.8 can be identified as ob¬ 
jects that drop-out in the i band. The constraint on the 
Y — J colour selects against cool, red stars that can be very 
faint in i. 

The selection criteria were applied to data obtained 
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Figure 1. The VIKING and KiDS filters used for the quasar 
search presented here plotted together with a (model) high- 
redshift quasar at 2 = 5.8, a quasar at 2 = 6.5 and one of 
the main contamin ants, a cool late-type star (data taken from 
iGeballe et Zll2n02h . Above a redshift of z ~ 5.8, the Lyman a 
line of the quasar shifts out of the i band, which results in a 
large measured i — Zv colour. Above approximately z ~ 6.5 the 
Lyman a line falls outside the Zy band filter. Although both 
high-redshift quasars and cool stars are very faint in the optical 
i band, the cool stars generally have a much redder Y — J colour 
as compared to quasars. 


in three areas NGP, SGP and GAMA09. The NGP and 
GAMA09 area are also covered by the SDSS. Where SDSS 
data were available, we applied additional colour criteria, 
based on the SDSS u (ua), SDSS g (ga) and SDSS i (is) 
catalogues. These additional criteria are: 

{is — Zy > 2.2 OR undetected in is) and undetected in its, 

9s, 

where undetected was defined as having a measured mag¬ 
nitude fainter than the 3cr limiting magnitude in that band 
(mag> mag3a). 

Due to the incompleteness of the KiDS i band cata¬ 
logues, especially at fainter magnitudes, and to objects close 
to the edge of the KiDS images, initially our selection criteria 
returned an unrealistically large number (several thousands) 
of potential z ^ 6 quasars. In the NGP and GAMA09 re¬ 
gion, the inclusion of the SDSS data significantly reduced 
the number of candidates, especially the number of candi¬ 
dates that were close to the edge of KiDS images. For re¬ 
maining candidates in those areas and those in the SGP area 
postage stamps in the i-band were made in order to reaffirm 
the magnitude of the sources. 

The postage stamps were analysed by measuring the 
flux of the candidates in an aperture centred on the VIKING 
position of the candidates. If a flux was measured with an 
S/N > 3, and the resulting colour was i — Zy < 2.2 the 
object was rejected as high-redshift quasar candidate. The 
remaining objects were visually inspected to confirm the re¬ 


ality of the candidate. A fraction was located at the edge, 
or just outside the KiDS imaging area and was removed 
from the list. We discovered several types of sources that 
mimic the colours of high-redshift quasars in our setup. Ob¬ 
jects moving between the epoch of the VIKING observations 
and KiDS observations and highly variable objects can have 
colours that are both blue in the near-infrared bands and 
very red in i — Zy. We were able to remove the majority of 
these objects as quasar candidates by carefully looking at the 
VIKING data, which are normally taken in different nights. 
Images in Zy, Y and J are generally taken in grey time, 
while images in J, H and Ka are taken during bright time. 
By comparing the two different J band epochs and rejecting 
objects with extremely blue Y — Ka colours {Y — Ka < —0.2, 
whi le quasars at 5.8 < a < 6.5 have 0.2 < Y — Ka < 0.5, 
e.g. iHewett et al.l[200^ ) we could remove most of the moving 
and variable objects from our list. 

In total, our final candidate list contained 30 sources: 6 
objects in the SGP region, 7 objects in the GAMA09 region 
and 17 candidates in the NGP region. 


3.2 Follow-up optical imaging 

We took images of 29 of our 30 2 ~ 6 quasar candi¬ 
dates during four separate observing runs. The first run was 
on 2012 June 29-30 at the 3.58 m ESO New Technology 
Telescope (NTT) using the European Southern Observa- 
tory’s Faint Object Spectrograph and Camera 2 (EFOSC2; 
iBuzzoni et al.l[l983l . We observed in the filters Gunn i (fil¬ 
ter i no. 705, hereafter 7n) and Gunn 2 : (filter 2 no. 623, 
hereafter 2 n). The second observing run was also with 
NTT/EFOSC2 on 2012 July 23-25. Thirdly, two c andidates 
were observed with ACAM llBenn, Dee fc Agocsl l2008f l on 
the William Herschel Telescope (WHT) on 2012 November 
10 in service mode, using the filters i and 2 . Finally, we used 
the NTT again on 2013 March 13-15 to follow-up candidates 
with EFOSC2. 

The candidates were observed with total exposure times 
between 720 and 2880 s. The integration time for each indi¬ 
vidual target depended on the sky transparency, seeing and 
expected brightness of the source. Data reduction steps in¬ 
cluded bias subtraction, flat fielding using sky flats and sky 
subtraction using the unregistered science frames. The im¬ 
ages were calibrated using bright, unsaturated point sources 
within the field of view. Since the ACAM i filter is very 
similar to the VST i band, the zero-point of the ACAM i 
images were directly derived from the KiDS magnitudes of 
point sources in the field. The zero-points of the NTT im¬ 
ages were computed using relations between 7 n and and 
SDSS a nd VIKING bands der ived from synthesized colours 
of stars (IVenemans et al.]l2013h . The results of the follow-up 
photometry are listed in Tabled 

The follow-up imaging rejected 26 of the 29 ob¬ 
served targets as potential 2 ~ 6 quasars. The ma¬ 
jority of rejected candidates have colours consistent 
with M dwarfs, see Fig. [5] The four objects that re¬ 
mained good quasar candidates were J083955.36-1-001554.2 
(hereafter J0839-f0015), J114833.18-7005642.3 (hereafter 
J1148-70056), J121516.88-7002324.7 (hereafter J1215-70023) 
and J032835.51-325322.8 (hereafter J0328-3253). Note that 
for J0839-I-0015 we did not obtain any follow-up photometry 
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Table 2. Followup photometry of 29 candidate 2 ~ 6 quasars 
in the KiDS+VIKING catalogues. All magnitudes are in AB and 
limiting magnitudes are 3 ct limits. 


Object name 

Filter 

Magnitude 

Telescope 

J020935.59-323206.5 

-^N 

>24.62 

NTT 


Zn 

>22.83 

NTT 

J032835.51-325322.8 

-^N 

21.14±0.02 

NTT 

J083635.31-000347.9 

i 

24.36±0.45 

WHT 


z 

22.51±0.21 

WHT 

J084810.50-1-021948.2 

i 

23.08±0.10 

WHT 


z 

22.94±0.28 

WHT 

J090647.34-004216.0 

In 

23.46±0.43 

NTT 

J091026.91-011137.2 

In 

22.24±0.03 

NTT 

J091131.58-001305.5 

In 

23.08±0.10 

NTT 

J092251.90-1-003245.3 

In 

22.75±0.05 

NTT 

J113549.91-I-015924.4 

In 

23.07±0.07 

NTT 

J114833.18-1-005642.3 

In 

22.67±0.05 

NTT 


Zn 

21.46±0.15 

NTT 

J115049.81-013830.2 

In 

23.20±0.06 

NTT 

J121144.38-I-005348.8 

In 

23.46±0.08 

NTT 

J121404.53-I-020220.2 

In 

22.38±0.10 

NTT 

J121516.88-I-002324.7 

In 

22.93±0.10 

NTT 

J122839.49-I-010253.7 

In 

22.13±0.02 

NTT 

J140552.33-I-030110.3 

In 

22.65±0.04 

NTT 

J141157.27-001613.5 

In 

22.92±0.12 

NTT 

J142054.44-I-021247.3 

In 

21 .88±0.02 

NTT 

J142300.93-I-011240.0 

In 

22.53±0.09 

NTT 

J142426.55-I-011304.1 

In 

22.75±0.04 

NTT 

J143428.30-010043.1 

In 

22.69±0.03 

NTT 

J144130.49-I-004119.7 

In 

22.52±0.03 

NTT 

J144922.61-004404.9 

In 

22.85±0.04 

NTT 

J145250.27-004047.8 

In 

23.22±0.05 

NTT 

J145442.58-020118.8 

In 

22.59±0.03 

NTT 

J224319.75-331506.8 

In 

>23.70 

NTT 


Zn 

>23.24 

NTT 

J225656.49-333910.1 

-^N 

22.10±0.05 

NTT 

J232355.24-303801.7 

-^N 

22.82±0.05 

NTT 


Zn 

22.07± 0.12 

NTT 

J234331.78-313627.5 

-^N 

21.30± 0.03 

NTT 


(see Section [4.211 . For the four remaining objects we took an 
optical spectrum to confirm their nature. 

3.3 Spectroscopic observations 

The spectroscopic observations of our four remaining quasar 
candidates were carried out with the FOcal Reducer/low 
dispersion Spectrograph 2 ('F0RS2: lAppenzeller et aDll998l i 
on the 8.2 m Very Large Telescope (VLT) Antu. For all 
observations we used a long slit with a width of l.Sarcsec 
and the grism ‘600z’ with a central wavelength of 9010 A 
and a resolution oi R = 1390. J0328-3253 was observed on 
2012 October 7 for 1782 s in clear conditions with a seeing 
around 0.8arcsec. On 2012 December 12 J0839+0015 was 
observed for 1782 s with a seeing ranging between 1.2arcsec 
and 1.7 arcsec. J1215+0023 was observed on 2013 February 
7 and 2013 March 5 for a total integration time of 5280 s. 
The seeing during the February observations varied between 

1.4 arcsec and 1.7 arcsec, while in March the conditions were 
better with a seeing around 0.7 arcsec. Finally, a spectrum 
of J1148+0056 was taken on 2013 April 13 for 2697 s in a 
seeing of 0.9 arcsec. All observations were setup in such a 



Z,-J [AB] 

Figure 2. I^ — Y versus Zy — J diagram illustrating the outcome 
of the follow-up NTT imaging of our 2 ~ 6 quasar candidates. The 
blue and green circles represent the colours of main sequence stars 
and M dwarfs, respectively. The yellow triangles are simulated 
colours of L dwarfs, and the connected black circles illustrate 
the colours of quasars at various redshifts 2 > 5.7. The crosses 
with error bars represent the candidates observed with the NTT. 
Based on the /n — Zy, ~ 7n ~ V and Zy — J colours most 

observed candidates were identified as foreground interlopers. The 
open squares are three objects that remained good candidates 
after the follow-up imaging and for which we obtained optical 
spectroscopy (see Section 13.31 1. Additional Zn photometry (see 
Table [21 not shown in this figure) confirmed that J1148+0056 
was a good high-redshift quasar candidate. Note that for a fourth 
good quasar candidate, J0839+0015, we did not obtain follow-up 
photometry (see Section US- 

way that a nearby bright star was included in the slit to 
provide a trace for the candidate spectrum. 

Daytime exposures of He, Ar and Ne lamps were ob¬ 
tained to provide the wavelength calibration. By measuring 
the wavelength of sky emission lines in our science spectra 
we determined that the wavelength calibration had a typical 
uncertainties between 0.2 and 0.5 A. For the flux calibration 
we used observations with a 5 arcsec long slit (taken the same 
night as the science observations) of one of the spectropho- 
tometric standard st ars Feige 110, HD49 798, L TT4816 and 
LTT7379 (IOkelll990l: rHamuv et al1ll992l. Il994h . Finally the 
flux calibrated spectra were scaled to the VIKING Zy band 
magnitudes to account for slit losses. 

The reduced spectra of the four quasar candidates are 
shown in Fig. |3 


4 FOUR QUASARS AROUND A REDSHIFT 
OF 6 

The four targets of our follow-up spectroscopy all show 
a blue continuum with a sharp break between 8300 and 
8500 A, characteristic of quasars at 2 ~ 6. To measure the 
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Figure 3. Discovery spectra from VLT/FORS2 of three z ~ 6 quasars in the KiDS+VIKING survey. The spectra of J0839+0015 and 
J1215+0023 are boxcar averaged over five pixels. For comparison, a spectrum of the sky is plott ed at the top. The positio n of several 
emission lines are marked. The spectra are overplotted with the composite quasar spectrum from IVanden Berk et af] I^OdJ). 


redshift of the newly identified quasars, we measured the 
redshifted position of emission lines in the spectra. When 
visible, the emission lines we used for this purpose are OI + 
Sill at Arest = 1305.42 A, Cii at Arest = 1336.6 A and SiIV 
+ O iv] at Arest = 1398.33 A. The central wavelength of the 
e mission lines are taken from the composite quasar spectrum 
bv IVanden Berk et al.l 1 I 2 OOIII . To compute the absolute mag¬ 
nitudes of the quasars at a rest-frame wavelength of 1450 A 
(Mi45o), we fitted a power-law slope in regions without emis¬ 
sion lines. The regions used were 1270 < Arest /A< 1295, 


1315 < Arest/A< 1325 and 1345 < Arest/A< 1370. Below we 
describe the new quasars in more detail. A summary of the 
properties of the quasars is given in Table [3] 


4.1 J0328-3253 at z = 5.86 

This quasar was selected as a high signal-to-noise ratio can¬ 
didate in our KiDS-|-VIKING catalogues. This source has 
a Zv band magnitude of Zv = 19.8 and blue infrared 
colours of Zv — V = —0.1 and Y — J = —0.1 (Table [3]). 
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The source was detected in the i band with a magnitude 
of i = 22.14 ± 0.06 and thus had a i — Zv = 2.3. Follow-up 
photometry with the NTT confirmed a break in the spectral 
energy distribution of the source shortwards of the Zv band. 
The In — Zy ~ 1-3, consistent with a quasar at z 5.8. 
VLT/FORS2 spectroscopy confirmed that the source is at 
2 = 5.86 ± 0.03 (Fig. [3]). We measured the redshift using 
the emission lines O I + Si II, C II and Si IV + O iv] and 
found a redshift around z = 5.86 for all three lines. Our 
power-law fit to the continuum gives an absolute magnitude 
of Mi 45 o = —26.60 ±0.04. After our discovery of this quasar 
in 2012 , the object was independently select ed as a bright 
quasa r candidate in the VST-ATLAS survey (ICarnall et al.l 

[2^. 


4.2 J0839±0015 at z = 5.84 

This object was initially selected with i—Zy > 2.7, Zy—Y = 
—0.1 and Y — J = 0.3 (Tabled). Analysis of the i band 
postage stamp revealed a faint source with a magnitude of 
i = 23.46T0.33, resulting in an i — Zy ~ 2.4. The two epochs 
of near-infrared data did not show strong variability and the 
detection of a source in the i band ruled out that this object 
was moving. We therefore decided to take a spectrum of 
this source. Despite poor observing conditions the resulting 
spectrum clearly shows a break at ~8300 A (Fig. [3]). If this 
break is at the wavelength of Lyman q, then the redshift 
would be around z 5.83. This redshift was confirmed by 
measuring the position of faint emission lines. The Ol ± 
Sill gave a redshift of z ~ 5.81, while the Si IV ± Oiv] line 
gave z ~ 5.87. We therefore estimate that the redshift of 
this quasar is z = 5.84 ± 0.04. The continuum fit gives an 
absolute magnitude of M 1450 = —25.36 ± 0.11. 


4.3 J1148±0056 at z = 5.84 

The source J1148±0056 was selected as a faint (Zy = 21.8) 
candidate with colours i — Zy > 2.3, Zy — Y = —0.2 and 
Y — J = 0.4. The i band postage stamp showed no source 
at the position of the Zy band object. Follow-up imaging 
with the NTT measured a colour of In — = 1.2 and — 

J = —0.3, consistent with a high redshift quasar. Optical 
spectroscopy detected a source with continuum longwards 
of A 8300 A and faint emission lines at 8930 A and at 
9115 A. We identified these lines with Oi ± Sin and Cii 
redshifted to z = 5.84 ± 0.03. From a power-law fit to the 
continuum we estimate that the absolute magnitude of this 
quasar is M 1450 = —24.46T0.il, which makes it the faintest 
quasar in our sample. 


4.4 J1215T0023 at z = 5.93 

J1215T0023 was a high priority quasar candidate in our lists 
with colours i — Zy > 2.7, Zy — Y = —0.2 and Y — J = 
0.3. NTT imaging confirmed that the source had a break in 
the spectral energy distribution shortwards of the Zy band 
(7n — Zy = 1.5). From the FORS2 spectrum we measure a 
redshift of z = 5.93±0.03, based on the Ol ± Sill and SiIV 
± O iv] lines. From the continuum we measure an absolute 
magnitude of M 1450 = —24.67 ± 0.14. 



Mi450 [AB] 

Figure 4. Distribution of quasars in the redshift range 5.8 < z < 
6.4 as function of absolute magnitude. Plotted are the histograms 
of the 60 previously published quasars in this redshift range 
(white) and new quasars discovered with the KiDS/VlKlNG sur¬ 
vey pilot (green). The solid black curve indicat es the luminosity 
functi on for quasars at redshift 6 as derived bv iKashikawa et alj 
ll2ni.5h . using their ‘case 1’. The dashed green line indicates the 
selection limit, i.e., the intrinsically faintest quasars that we 
could select with the KiDS/VIKING survey. The aim of the 
KiDS/VIKING survey is to triple the number of known quasars 
fainter than M 1450 > —25. 


5 DISCUSSION AND OUTLOOK 


By combining KiDS i band data with VIKING 
Zy,Y,J,H,Ks catalogues, we performed a pilot study 
to search for faint z ~ 6 quasar candidates. In our first 
year of obtaining follow-up data, we discovered four 
new quasars in an area of ~254 deg^. To compare our 
quasars with those of other high redshift quasar searches 
we pl ot in Fig. |4] the z ~ 6 quasar luminosity function 
from iKashikawa et all ll2015h as function of absolute 
magnitude. Overplotted is a histogram of the 60 known 
quasars at 5.8 < z < 6.4, taken from t he literature 
|Zeimann_et_^lJ_ 201 l|; Banados_et_^ 2014 : ICalura et al ' 


2 OI 4 I: IKashikawa et al l[2015l:ICarnall et al.ll2015l:l,Iiang et al 


2 OI 5 I : Reed et aP l2015l . and references therein). The new 
quasars from the KiDS/VIKING surveys are generally 
fainter {Zy mag > 21) than the quasars found in, for 
example, the SDSS and PSl surveys, showing that our 
z ~ 6 quasar search is probing the faint end of the known 
z = 6 quasar population. 

One of the questions is how many quasars we were ex¬ 
pecting in our survey. To estimate this, we took the z ~ 6 
quasar luminosity function and computed the expected num¬ 
ber of quasars as function of area and limiting magnitude. 
The limiting magnitude varied from area to area, and mostly 
depended on the depth of the KiDS i-band imaging. The 
median depth to which we coul d search for q uasars was 
Zy ~ 21.65. We integrated the IWillott et al.l ll2010bl ) lu¬ 
minosity functior0 and in the redshift range 5.8 < z < 6.4 


® Using_tliejKashikaw^_etalJ ||2015^ luminosity function instead 
of the I Willott ^t’afTI^niObh luminosity function gives the same 
result. 
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Table 3. Properties of the newly discovered quasars. The redshift uncertainty does not include the uncertainty between the redshift 
of the UV emission lines and the systemic redshift of the quasars. The uncertainty in the absolute magnitude is a combinations of the 
uncertainty in the continuum fit and in the absolute scaling of the spectrum based on the Zy-band magnitude. 


Object name 

RA (J2000) 

Dec. (J2000) 

Redshift 

•^VIKING 

UviKING 

•IviKING 

Ks,VIKING 

A7i450 

J0328-3253 

J0839-I-0015 

J1148-1-0056 
J1215-I-0023 

03''28™35!511 
08*139™ 55? 356 
ll‘'48™33?180 
12‘'15™16!879 

-32°53'22792 

+00‘’15'54'.'21 

+00‘’56'42'.'26 

+00‘’23'24'.'74 

5.86±0.03 

5.84±0.04 

5.84±0.03 

5.93±0.03 

19.83±0.02 

21.09±0.05 

21.79±0.10 

21.42±0.07 

19.89±0.04 

21.19±0.10 

21.96±0.23 

21.60±0.18 

19.98±0.03 

20.89±0.11 

21.55±0.24 

21.32±0.14 

19.69±0.05 

20.50±0.15 

21.20±0.23 

21.10±0.21 

-26.60±0.04 

-25.36±0.11 

-24.46±0.11 

-24.67±0.14 


we expected ~7 quasars in the area we surveyed, where we 
discovered four. The difference between the predicted num¬ 
ber and the observed number of quasars might well be due 
to various sources of incompleteness of our search. The most 
prominent ones are: catalogue completeness, point source se- 
lection efhciency and c ompleteness of the colour selection. In 
IVenemans et alj ll2013ll it was found that the VIKING cata¬ 
logue completeness is ~96% above a S/N of 7, while the frac¬ 
tion of point sources that is incorrectly classihed as extended 
objects is around 6%-7%. The largest incompleteness is the 
colour incompleteness: both due to random photometric er¬ 
rors and due to small variations in the quasar spectral energy 
distribution and in the intergalactic medium (IGM) absorp¬ 
tion, quasars can be placed outside our colour selection box. 
For example, quasars with a redshift around 2 ~ 5.8 have an 
i —Zy colour that is a st eep function of redshift, see e.g. fig. 1 
in lBafiados et all (l2014h . It is therefore possible that we are 
missing quasars around 2 ~ 5.8 that have i — Zy < 2.2. In¬ 
deed, by matching our data with high-redshift quasars that 
have previously been discovere d, we found tha t the SDSS 
quasar J0836-I-0054 at 2 = 5.82 dFan et al.ll200lh is detected 
in both the VIKING and KiDS survey. While the redshift 
of this quasar is very close to the ones we present in this 
paper, J0836-I-0054 has an i — Zy = 1.70, well below our 
cut of i — Zy = 2.2. At the high end of our probed red¬ 
shift range, 2 ~ 6.4, due to variations in the IGM and in 
the quasar’s emission line strength and due to photomet¬ 
ric uncertainties, quasars can become too red in Zy — V to 
be selected b y our colour crit e rion ( Zy — Y< 1.0, see Sec¬ 
tion [3T]). In Venemans et al.l (l2013ll it was shown through 
modelling that up to 15% of artihcial quasars at 2 = 6.4 
fulfilled the colour selection Zy — Y > 1.1. This indicates 
that we are also incomplete in our selection of quasars at 
2 = 6.3 — 6.4. Gombined with the other sources of incom¬ 
pleteness, the discovery of four new quasars is well in line 
with the theoretical expectation of seven. 


The total area that ultimately will be covered by the 
KiDS and VIKING surveys is 1500 deg^, roughly six times 
the area we searched in our pilot study. We therefore ex¬ 
pect to discover at least 20 additional quasars in this area. 
Furthermore, there are several ways to push our quasar 
search to even fainter magnitudes. For example, when se¬ 
lecting quasar candidates, we could reduce the i — Zy colour 
lim it to i — Zy = 2.0 (similar to the PSl quasar search, 
see iBanados et al.l [2oT3) . This will increase our complete¬ 
ness, especially at the low-redshift end of our quasar search, 
2 5.8. Additionally, our requirement for a candidate to 

have a detection in the Y band, while significantly reducing 
the number of spurious sources, limits in certain areas how 
faint we can select candidates as the i and Zy band cata¬ 


logues are deeper than the V-band catalogue (see Table [T]). 
By pushing our selection to fainter magnitudes, we should 
be able to uncover quasars as faint as M 1450 = —24. Based 
on the quasar luminosity function, we therefore expect to 
discover at least 30 more quasars down to M 1450 = —24, 
significantly increasing the number of quasars at these faint 
absolute magnitudes. 

Besides exploring the new survey data that will become 
available, our plan is to characterize the newly discovered 
quasars with near-infrared, far-infrared and sub-mm spec¬ 
troscopy. This will allow us to constrain black hole masses, 
ISM metallicity and host masses. This is important to un¬ 
derstand the relationship between the accreting black hole 
and the stellar component of galaxies in the first billion years 
of the Universe. 
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